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Abstract: In this paper we discuss the effect of the complete leading-order QED correc¬ 
tions to the DGLAP equations that govern the perturbative evolution of parton distribution 
functions (PDFs). This requires the extension of the purely QCD DGLAP evolution includ¬ 
ing a PDF for the photons and, consistently, also for the charged leptons e^, and r^. 
We present the implementation of the QED-corrected DGLAP evolution in the presence of 
photon and lepton PDFs in the APFEL program and, by means of different assumptions for 
the initial scale PDFs, we produce for the hrst time PDF sets containing charged lepton 
distributions. We also present phenomenological studies that aim to assess the impact of 
the presence of lepton PDFs in the proton for some relevant SM (and BSM) processes at 
the LHG at 13 TeV and the FCC-hh at 100 TeV. The impact of the photon PDF is also 
outlined for those processes. 

Keywords: parton distribution functions, electroweak corrections, high-precision compu¬ 
tation 
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1 Introduction 

The Run II at the Large Hadron Collider (LHC) will probe the interactions of elementary 
particles at unprecedented scales of energy and with an outstanding accuracy. In order 
to correctly identify possible Beyond-the-Standard-Model (BSM) effects, the theoretical 
predictions for the Standard Model (SM) processes have to match (at least) the precision 
of the corresponding experimental measurements. In other words, the impact of higher- 
order corrections on phenomenological predictions has to be known and under control. 
To this purpose, the computation of Next-to-Leading-Order (NLO) QCD corrections is 
necessary, but often not sufficient. For example, at fixed order the inclusion of the second- 
order (NNLO) corrections in QCD as well as of the NLO electroweak (EW) corrections 
is in general desirable and in particular cases even essential. In order to formally achieve 
this level of accuracy, not only the matrix elements of the hard processes, but also the 
so-called parton-distribution functions (PDFs) have to be known at the same level of pre¬ 
cision. While most of the collaborations already provide PDF sets accurate up to NNLO 
in QCD [1-5], for the EW corrections the situation is less satisfactory. Indeed, NLO EW 
corrections naturally introduce photon-induced processes and only two PDF sets includ¬ 
ing a photon density are presently available: the MRST2004QED set [6], whose photon 
PDF is determined by means of model assumptions, and the NNPDF2.3QED family [7], 
where instead the photon PDF is extracted by means of a fit to data. Contrary to the 
MRST2004QED set, the photon PDF of NNPDF2.3QED sets is provided with an uncer¬ 
tainty. In addition, the aforementioned PDF sets are accurate up to (N)NLO in QCD but 
only at LO in QED evolution. Thus, since NLO EW, and in particular NLO QED effects. 
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are never included at the level of PDFs, this level of accuracy cannot be formally claimed 
at the level of the hadronic cross section. 

Before moving to NLO QED accuracy, however, it is worth attempting a determination 
of PDFs at (N)NLO QCD + LO QED accuracy for all the QED-interacting SM particles, 
i.e. by considering also leptons as partons of the proton. In this way predictions for lepton- 
induced processes at hadron colliders could be easily calculated by convoluting partonic 
cross sections with lepton PDEs. On the one hand, this can be useful for the study 
of processes such as same-sign and/or different-flavour dilepton production, which, with 
lepton PDFs, can be directly produced without requiring additional (unresolved) jets or 
photons. On the other, lepton-initiated processes are naturally involved in EW corrections. 
Indeed, as already mentioned, a very interesting aspect of EW corrections is the fact that 
they naturally require new partons in the proton on top of the usual (anti)quarks and 
gluons. Processes with quarks in the initial state at LO receive NLO EW corrections 
from photon-initiated processes. Equivalently, once the photon is considered as a parton, 
photon-initiated processes at LO will receive corrections from processes including leptons 
in the initial state. However, while PDF fits that provide an estimate of the photon PDFs 
are presently available, this is not the case for leptons. 

There are (at least) two reasons why lepton PDFs have been neglected so far in the 
context of EW corrections. Eirstly, they are expected to be small, since they are suppressed 
by a factor a with respect to the photon PDE, which in turn is already suppressed with 
respect to quark and gluon PDEs. Secondly, photon-initiated processes are typically sim¬ 
ulated only at LO accuracy, due the smallness of the photon PDF with respect to quark 
and gluon ones. Thus, neglecting NLO EW corrections in photon-initiated processes, the 
absence of the lepton PDEs has never been a real issue. On the other hand, the recent 
progresses in the automation of NLO EW corrections [8-12] make it possible to take into 
account these contributions without additional efforts. Moreover, lepton PDFs can induce 
not only EW corrections, but also new LO contributions for many SM and BSM processes 
at hadron colliders. In particular, leptons in the initial state can open new production 
mechanisms. Eor instance, at LO the production of a lepton pair can proceed also through 
a t-channel photon exchange as it is the case for the elastic scattering at lepton colliders. 
This production mechanism is enhanced in the peripheral region close to the beam pipe, 
therefore a quantitative estimate of lepton PDFs is desirable in order not to introduce any 
systematic effect in the searches and measurements in the Run II at the LHC and possibly 
also at future colliders. 

The purpose of this paper is to give the first quantitative estimate on the leptonic 
content of the proton and analyse its phenomenological impact at the LHC at 13 TeV 
and at a future 100 TeV hadron-hadron collider (FCC-hh). This will be achieved in two 
steps. First of all by implementing the lepton PDF Dokshitzer-Gribov-Altarelli-Lipatov- 
Parisi (DGLAP) evolution equations at LO in QED in the so-called variable-flavour-number 
(VEN) scheme in the evolution code APFEL [13], and secondly by producing a guess for the 
lepton PDFs at the initial scale based on the assumption that leptons are generated by 
photon splitting at the respective mass scales. 

The paper is organised as follows. In Section 2 we describe the procedure used to 
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generate new PDF sets including leptons. First, in Section 2.1 we show how to generalise 
the DGLAP evolution equations in order to include the complete QED corrections at LO, 
which involve leptons. Then, in Section 2.2 we discuss how to model the lepton distributions 
at the initial scale on the base of different theoretical assumptions. In Section 2.3 we 
present the numerical results obtained with APFEL, considering different initial conditions 
and comparing the results at the level of PDFs and momentum fractions. In Section 3 we 
discuss the phenomenological impact of lepton (and photon) PDFs at the LHC and at the 
FCC-hh. In Section 3.1 we look at the impact of lepton PDFs on parton luminosities and 
we analyse their dependence on the invariant mass and the rapidity of the initial/final state 
in a process-independent approach. Then, in Section 3.2, we turn to the study of the effects 
of lepton PDFs for some relevant SM (and BSM) processes at the LHC and at the FCC- 
hh. Finally, in Section 4, we draw our conclusions and discuss possible phenomenological 
studies that could be done using PDF sets with lepton distributions. The PDF grids used 
in this work have been made public in the LHAPDF6 library [14] format. 


2 PDFs determination 


2.1 DGLAP equations in the presence of photons and leptons 


In this section we will show how to extend the DCLAP equations in order to include the 
evolution of photon and lepton PDFs at leading order (LO) in QED. The inclusion of the 
photon PDF evolution has already been treated extensively in Refs. [ 6 , 13, 15] and thus 
we will not discuss it here. As far as leptons are concerned, at LO in QED they couple 
directly only to photons. However, since the photons couple to quarks that in turn couple 
to gluons, the lepton PDFs evolution will indirectly depend on the evolution of all other 
partons. Following the notation of Ref. [13], where QCD and QED evolutions are treated 
separately, the inclusion of leptons does not imply any change to the pure QCD sector. On 
the contrary, with the inclusion of leptons, the QED evolution equations with respect to 
the QED factorisation scale v read: 
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where 7 , q{q) and £ (£+) correspond to the PDFs of the photon, the quark(antiquark) 
flavour q and the lepton(antilepton) species respectively. In Eq. (2.1), corresponds to 
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the electric charge of the quark flavour q, Nc = 3 is the number of QCD colours and o. 
is the running fine-structure constant. Note also that the indices q and I in the sums in 
the first line of Eq. (2.1) run over the Uf and active quarks and leptons at the scale u, 
respectively. The symbol (8> represents the usual Mellin convolution operator defined as: 

dy f dz A{y)B{z)5{x — yz). (2.2) 

Jo 


A{x) ® B{x) = 


Finally, the LO QED splitting functions are given by: 

= -^5(1-x), 

Pq^Kx) = P^!^\x) = 2[x^ + {1- x)2] , 
P\yx) = P\^^{x) = 2\^-±^^^^ , 

Pm^^) = Pu\^) = ^ (1 ^x)_^ + • 


(2.3) 


Combining the system of differential equations in Eq. (2.1) with the pure-QCD DGLAP 
equations that govern the evolution of gluons and quarks, we obtain the full (N)NLO QCD 
plus LO QED evolution in the presence of photons and leptons. The solution of this system 
is implemented in version 2.4.0, or later, of the APFEL evolution library [13]. 

It is worth stressing that, contrary to the approach described in the original reference 
of APFEL [13], in the new implementation of the coupled QCD-I-QED DGLAP equations 
the QCD and the QED evolutions are no longer treated separately and combined a poste¬ 
riori. Indeed, thanks to a suitable combination of PDFs that cleverly diagonalises the full 
splitting-function matrix in the presence of QED corrections, we have been able to imple¬ 
ment in APFEL an efficient unified solution of the QCD-I-QED DGLAP equations where the 
evolution of PDFs with respect to QCD and QED factorisation scales takes place at the 
same time. In other words, QCD and QED factorisation scales are no longer kept separate 
and are instead identified. We present the details of this basis in Appendix A. 

As far as the fine-structure constant a is concerned, consistently with the evolution of 
PDFs, we only consider the LO running by solving the renormalisation-group equation: 

= /5qed«^(j^), (2.4) 

with: 

where, as a boundary condition for the evolution, we take a~^(mT-) = 133.4. 
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2.2 Modelling the lepton PDFs 


Once the full QCD+QED DGLAP evolution in the presence of photons and leptons has 
been implemented, the following step is the determination of suitable boundary conditions, 
i.e. the initial scale PDFs, to be evolved. As far as gluons and quarks are concerned, 
various solid determinations are publicly available [1-5]. On the other hand, as already 
mentioned, only two public determinations of the photon PDF are available at the moment: 
the MRST2004QFD set [6], where the photon PDF is determined by means of a model 
assumption, and the more recent NNPDF2.3QFD family [7], where the photon PDF is in¬ 
stead determined from a fit to data. No lepton PDF determination has ever been attempted 
so far. 

In principle, the most reliable way to determine the lepton PDFs is the extraction from 
data by means of a fit. However, to our opinion, this option is currently not viable. The 
main reason is that lepton PDFs are expected to be tiny as compared to the quark and gluon 
ones and even substantially smaller than the photon PDF. Assuming no intrinsic leptonic 
component in the proton, lepton PDFs are expected to be of the order of a times the 
photon PDF, where a ~ 10“^. As a consequence, the photon being already suppressed as 
compared to quark and gluon PDFs, the contribution of leptons is likely to be so small that 
the determination of the corresponding PDFs from experimental data is hardly achievable. 

As an alternative to the fit, one could try to guess the functional form of the lepton 
PDFs just by assuming that light leptons, i.e. electrons and muons, are generated by photon 
splitting at the respective mass scales. Using the leading-logarithmic (LL) approximation 
for the photon splitting into leptons, one can then infer their distributions at the initial 
scale Qo — 1 GeV as: 


e {x,Qo) = £'^{x,Qo) 


47r \fnj J 


^p(O) 


l{y,Qo), 


( 2 . 6 ) 


with £ = e, n ^ . In other words, Fq. (2.6) comes from the linearisation of the LO QFD 
DGLAP evolution for the lepton PDFs from the mass scale m^, where they are assumed 
to appear as dynamical partons, to the initial scale Qq. Indeed, in the situation where 
both QFD and QGD evolutions are considered, the divergence of the strong coupling as in 
correspondence of Aqcd (Landau pole), with Aqcd > makes such evolution impossible 
to implement in an exact way and Eq. (2.6) represents an approximated solution. Thus, 
Eq. (2.6) is equivalent to assume that the photon PDF does not evolve between the lepton 
mass scale rrii and the initial scale Qq. However, the DGLAP evolution between m£ and 
Qo would introduce additional 0{a) terms that exceed the LO QFD accuracy. 

Since the tau mass is typically larger than Qq {rrir = 1.777 GeV > Qq), we have chosen 
to determine the PDFs according to the usual VFN scheme [17-19], i.e. by dynamically 
generating them at the threshold. 
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ID 

PDF Set 

Ref. 

QCD 

QED 

Photon PDF 

Lepton PDFs 

A1 

apfel nn23nlo0118 lept0 

[20] 

NLO 

LO 

l{x,Qo) = 0 

Eq. (2.7) 

A2 

apfel nn23nnlo0118 lept0 

[20] 

NNLO 

LO 

O 

II 

o 

Eq. (2.7) 

B1 

apfeljnn23qedlo0118 lept0 

[21] 

LO 

LO 

Internal 

Eq. (2.7) 

B2 

apfel nn23qednlo0118 lept0 

[7] 

NLO 

LO 

Internal 

Eq. (2.7) 

B3 

apfel nn23qednnlo0118 lept0 

[7] 

NNLO 

LO 

Internal 

Eq. (2.7) 

B4 

apfel mrst04qed lept0 

[6] 

NLO 

LO 

Internal 

Eq. (2.7) 

Cl 

apfeljnn23qedlo0118 lept 

[21] 

LO 

LO 

Internal 

Eq. (2.6) 

C2 

apfel nn23qednlo0118 lept 

[7] 

NLO 

LO 

Internal 

Eq. (2.6) 

C3 

apfel nn23qednnlo0118 lept 

[7] 

NNLO 

LO 

Internal 

Eq. (2.6) 

C4 

apfel_mrst04qed_lept 

[6] 

NLO 

LO 

Internal 

Eq. (2.6) 


Table 1: Summary of the sets of PDFs including photons and leptons generated with APFEL. 


2.3 PDF sets with lepton distributions 

In this section we discuss the results of the implementation of the lepton PDF evolution 
in APFEL. As discussed in the previous section, the determination of lepton PDFs from a 
direct fit to data seems to be hard to achieve and thus, as an alternative, we adopt a model 
based on a theoretical assumption for modelling lepton PDFs at the initial scale. 

The model presented in the previous section is based on the assumption that leptons 
are generated in pairs from LL photon splitting at the respective mass scales. This results 
in the ansatz in Eq. (2.6) for the light lepton PDFs. However, in order to test how sensitive 
the results are with respect to the initial scale distributions, we also consider the “zero- 
lepton” ansatz where the lepton PDFs at the initial scale Qq are (artificially) set to zero, 
that is: 

r{x,Qo) = £^{x,Qo) = 0. (2.7) 

In this context, the construction of PDF sets with leptons requires a pre-existing PDF 
set to which we add our models for the lepton distributions. Of course, in order to apply 
the ansatz in Eq. (2.6), we need PDF sets that already contain a photon PDF, e.g. the 
MRST2004QED set [6] or the NNPDF2.3QED sets [7]. We will use both sets to generate 
lepton PDFs via the ansatz in Eq. (2.6). On the contrary, the ansatz in Eq. (2.7) can be 
applied to any set so that lepton and photon distributions are generated just by DGLAP 
evolution. 

In order to assess the effect of considering lepton PDFs in the DGLAP evolution, we 
consider three different initial-scale configurations that are also summarised in Tab. 1: 

• sets where both photon and lepton PDFs are initially absent and dynamically gen¬ 
erated by DGLAP evolution. For this configuration we have constructed the sets A1 
and A2 in Table 1 based on NNPDF2.3 NLO and NNLO, respectively. 

• Sets where the photon distribution is present and the lepton PDFs are set to zero at 
the initial scale Qq {i.e. Eq. (2.7)) and then evolved as discussed in Sect. 2.1. These 

^For the light lepton masses, we take rrie = 0.510998928 MeV and mi_^ = 105.6583715 MeV, as quoted 
in the PDG [16]. 
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configurations are based on the NNPDF2.3QED and MRST2004QED sets of PDEs 
and identified by the indices Bl, B2, B3 and B4 in Table 1. 

• Sets of PDEs generated from NNPDE2.3QED and MRST2004QED using the ansatz 
in Eq. (2.6) for the light lepton PDEs (sets Cl, C2, C3 and C4 in Table 1). 

The evolution of the PDE sets listed in Tab. 1 is performed using APFEL, as discussed 
in Sect. 2.1. All PDE sets are tabulated in the LHAPDF6 format, which allows us to include 
the lepton PDEs in a straightforward manner. In the following we will present the results of 
the implementation of the lepton PDE evolution for the different ansatze discussed above. 
In particular, we will quantify momentum fractions, correlations and the differences among 
the different sets. In order to make the discussion more concise, we will only presents results 
for the NLO sets. However, the features we observe are common to all sets and thus all 
the conclusions we draw apply also to LO and NNLO sets. 

In Eig. I we show the lepton and photon PDE central values for the set identified by 
the ID AI in Table 1. In this configuration photons and leptons are set to zero at Qo = 1 
GeV and then dynamically generated by DGLAP evolution. The left plot shows PDEs at 
Q = 1.8 GeV, in this case electron and muon PDEs are identical (by definition), and the r 
PDE has just been dynamically generated {rrir = 1.777 GeV). In the right plot the same 
comparison is displayed at Q = 100 GeV, showing that all lepton PDEs are close to each 
other. Similar results are obtained also with the NNPDE2.3 NNLO (A2) set. 

Configurations B2 and B4 are shown in Eig. 2. In this case, the prior sets of PDEs con¬ 
tain a photon PDE while lepton PDEs are null at the initial scale Qo and thus generated dy¬ 
namically by DGLAP evolution. Again, similar results are obtained for the NNPDE2.3QED 
LO (Bl) and NNLO (B3) sets. 

Einally, we consider the configurations of type C where, starting from a prior containing 
a photon PDE, the initial distributions for the leptons are determined using the ansatz in 
Eq. (2.6). In Eig. 3 we show the resulting lepton PDEs for the configurations C2 (top) 
and C4 (bottom), at Q = 1.8 GeV (left) and Q = 100 GeV (right). Again, the qualitative 
behaviour for the configurations Cl and C3 is the same. 

Eig. 4 shows the ratios for the light lepton PDEs at Q = 100 GeV for the configurations 
Al, B2, C2 to G2. In this way we can quantify the difference generated by the various 
initial conditions starting from the same prior set, i.e. NNPDE2.3 at NLO. Eor the electron 
PDEs (left plot), the ansatze in Eqs. (2.6) and (2.7) applied to a set with a photon PDE 
leads to similar results in the small-x region, while differences up to 50% are observed in 
the large-x region. The electron PDEs resulting from the set without a photon PDE are 
instead much smaller for all in the entire x range. The same behaviour is observed also for 
the muon PDEs (right plot in Eig. 4), with slightly less enhanced discrepancies as compared 
to electrons. 

As already mentioned, the photon PDE of the NNPDE2.3QED sets is provided with 
an uncertainty. This enables us to assess how compatible the lepton PDEs emerging from 
the three different configurations described above are. In Eig. 5 we show the same plots 
presented in Eig. 4 but including the 68% confidence level (CL) bands (when available). We 
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Figure 1: Lepton and photon PDFs for the configuration Al at Q = 1.8 GeV (left) and 
Q = 100 GeV (right). 






Figure 2: Same as Fig. 1 for the configurations B2 (top) and B4 (bottom). 


observe that the configurations B2 and G2 are perfectly compatible within uncertainties 
all over the range in x. We can thus conclude that, even though the ansatz in Eq. (2.6) 
provides an arguably better estimate of the lepton PDFs, the present precision with which 
the photon PDF is determined is such that the difference with respect to the “zero-lepton” 
ansatz in Eq. (2.7) is still within uncertainties. Finally, it is unsurprising that the con¬ 
figuration Al, that has no uncertainty band because it assumes that the photon PDF is 
identically zero at the initial scale for all PDF replicas, lies right on the lower limit of the 
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Figure 3: Same as Fig. 1 for the configurations C2 (top) and C4 (bottom). 




Figure 4: Electron (left) and muon (right) PDFs at Q = 100 GeV for the configurations 
Al, B2 and C2, displayed as ratios to C2. 


uncertainty band of both configurations B2 and C2. The reason for this can be traced 
back to the fact the the photon PDF of the NNPDF2.3QED sets at the initial scale is 
compatible with zero and, due to the enforcement of positivity, the lower limit of the 68% 
CL band essentially concides with 7 (x, Qq) = 0. 

Interesting information about the photon and lepton content of the proton is provided 
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Figure 5: Same as Fig. 4 but with 68 % CL bands displayed. 


by the respective momentum fractions defined as: 

MF^((5) = f dxx''y{x,Q), and MF£±((5) = f dxxi^{x,Q). (2-8) 

Jo Jo 

In Fig. 6 we plot the percent momentum fractions as a function of the scale Q for the 
conhgurations B2 (left) and C2 (right). While the photon PDF carries up to around 1% 
of the momentum fraction of the proton, lepton PDFs carry a much smaller fraction (by 
about two orders of magnitude), regardless of the initial-scale parameterisation. This is 
consistent with the fact that, for both parameterisations, lepton PDFs are proportional 
to a times the photon PDF (£ oc a x 7 ). In conclusion, lepton PDFs carry such a small 
fraction of the proton momentum that they do not cause a signihcant violation of the total 
momentum sum rule when attached to a pre-existing PDF set. 

Now we turn to consider the uncertainties of the lepton PDFs. A realistic estimate 
of the lepton PDFs also requires an estimate of the respective uncertainties. To this end, 
using exactly the same procedure discussed in the previous sections, we have generated 
lepton PDFs for all replicas of the NNPDF2.3 family sets. This allows one to estimate the 
uncertainty on each lepton PDF in the usual way. In Fig. 7 we plot the lepton PDFs with 
the respective uncertainty for the configurations C2 and C3 at Q = 100 GeV. Uncertainties 
are presented as symmetric 68 % CL bands centred around the mean value of each PDF. As 
expected, the lepton PDF uncertainties follow the pattern of the photon PDF, characterised 
by a large uncertainty at large x. 

To conclude this section, in Fig. 8 we present the PDF correlation matrix for the 
conhguration C3 at Q = 1.8 GeV on a grid of = 50 points logarithmically distributed 
in the interval xi,X 2 = [10“®, 0.9], for the partons r, jj., e, 7 , g, d, u and s. The correlation 
coefficients are defined as: 




Wep / Q)fl3{x2, Q))^gp - {fajxi, Q)) 

Arep - 1 I Cra{xi,Q) ■ ai3{x2,Q) 


(2.9) 

where averages are taken over the Aj-ep replicas and where ai{x,Q) are the corresponding 
standard deviations. 
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Figure 6: Momentum fractions for the photon and lepton PDFs. 




Figure 7: Uncertainty bands for the photon and the lepton PDFs of configurations C2 
(left) and C3 (right) at Q = 100 GeV. 

We note a clear distinction between the QED (top-left square region) and the QCD 
sector (bottom-right region). As expected, there are strong correlations between r, /x, e 
and 7 due to the fact that leptons are generated by photon splitting. A similar but milder 
behaviour is also observed for g, d, u and s. Finally, the off-diagonal elements show that 
quark and gluon distributions are instead very mildly correlated to lepton and photon 
PDFs. Similar results are obtained also for the other configurations. 

3 Phenomenological impact 
3.1 Parton luminosities 

In the computation of the cross sections of hadron-collider processes, PDF contributions 
factorise in the form of parton luminosities. Thus, before looking at specific processes, it 
is useful to study the behaviour of the parton luminosities of the different initial states, by 
including also leptons. 

Parton luminosities are doubly differential quantities defined as: 
fpr.. 

= fiixi,Q) fj ix 2 ,Q) , Xl = ^/Tey, X2^VTe~y, T = Mxjs, (3.1) 
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Figure 8 : PDF correlation coefficients at Q = 1.8 GeV for the configuration C2. 


where s is the squared center-of-mass energy of the hadronic collision, and Mx and y 
are the invariant mass and the rapidity of the partonic initial/final state, respectively. 
In Eq. (3.1), fi{x,Q) is the PDF of the r-th parton evaluated at the factorisation scale 
Q. Different choices for Q can be adopted in order to improve predictions of a particular 
process and/or distribution. At the level of pure luminosities, without the convolution with 
any specific matrix element, the factorisation scale can be naturally set to Q = Mx- 

For the numerical evaluation of the different luminosities, it is practical to integrate 
out one of the kinematic variables in Eq. (3.1), namely, either r or y. In particular, by 
integrating out y, one obtains the Mx-differential luminosities, which read: 

dC ■ ■ 1 cIt 

{Mx) = = - —ft (x, Mx) fj {tIx, Mx) - (3.2) 

By integrating out r, instead, the y-differential luminosities are given by: 
dC- 

'^ijiy) = = 2e“^^ / dxxfi{x,^/sxe~'^)fj{xe~‘^y,^/sxe~y ), (3.3) 

J yjTent eV 

with Tcut = ^Xenth- (3-3) the lower bound of the integral, proportional to y/rfff, 

implies that Mx > Mx,cut- Although artificial, this kind of cut is always present in 
phenomenological predictions. Indeed, in the computation of hadronic cross sections a cut 
on Mx is automatically introduced by the masses of the particles in the final state or, 
indirectly, by the cuts on the massless particles, which are necessary to avoid infrared and 
collinear divergencies and ensure the finiteness of the cross sections. 

In order to compare the size and the Mx-dependence of the different parton luminosi¬ 
ties, in Fig. 9 we plot ^iq, ^qg, ^gq and ^gg as functions of 
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Figure 9: Mx-differential parton luminosities involving also photon and lepton PDFs. 


Mx at y/s = 13 TeV (left) and y/s = 100 TeV (right) for the set C2. The quantities 
<!)£+£-, ^iq and are dehned as: 

{Mx) = {Mx) , {Mx) = {Mx) , (3.4) 

j=e=‘=,/i±,T± j=e,fi,T 

^iq{Mx)^ Y, ^ijiMx), %g{Mx)^ Y ^ni^x). (3.5) 

j=u,u,...,b,b j=u,d,c,s,b 

Similarly, in Fig. 10 we plot the Ty.g, Ty.g, Tgg, Tgg and Tgg parton 

luminosities as functions of y at the same energies and for the same PDF set. In the case 
of y/s = 13 TeV (left), we set Mx,cut = 10 GeV and for y/s = 100 TeV (right) we set 
Mx,cut = 100 GeV. The Tjg and Tgg luminosities are defined in full analogy 

with Eqs. (3.4) and (3.5). In Figs. 9 and 10 we also plot with consistent colours the 68% 
uncertainty of the luminosities involving photon or lepton PDFs (lower insets) and of the 
luminosities involving only (anti)quarks or gluons (central insets). 

The relative size of the plotted luminosities follows the expected pattern. In general, 
the photon PDF suppresses the luminosity by a factor of a w.r.t. the (anti)quark PDFs 
and, analogously, the lepton PDFs suppress the luminosity by an additional factor of o. 
w.r.t. the photon PDF. This can be easily seen in Fig. 9(10), e.g. by comparing <^^,£(4'^^) 
with 4>y,y,(4'^y,) und (4'£+£-), the three lowest curves. However, from Fig. 9 we also no¬ 
tice that this hierarchy is not clearly respected at large invariant masses. In this kinematic 


- 13 - 
























apfel_nn23qednlo0118_lept (C2) apfel_nn23qednlo0118_lept (C2) 



y y 


Figure 10: y-differential parton luminosities involving also photon and lepton PDFs. 


region the QCD luminosity combinations, and get closer to the luminosities 

involving photons and leptons, suggesting possibly non-negligible phenomenological impli¬ 
cations due to lepton and photon channels. The major part of this effect is caused by the 
relative behaviour of the strong coupling as with respect to the QED coupling a as func¬ 
tions of the scale Mx- As is well known, DGLAP evolution leads to an increase of PDFs in 
the small-x region and to a decrease in the large-x region as the evolution scale increases. 
The magnitude of such a behavior is driven by the rate of change of the coupling, no matter 
whether it is negative, as for as, or positive, as for a. Given that the absolute value of the 
QGD /3-function is larger than the QED one up to very large scales, a^-driven PDFs, like 
quark and gluon PDFs, are relatively more suppressed at larger values of x as compared to 
a-driven PDFs, like lepton and photon PDFs. It is clear from Eq. (3.2) that the behavior 
of the Mx-differential luminosities at large values of Mx reflects the behavior of PDFs 
at large values of x and thus, following the argument given above, QGD luminosities are 
more suppressed than QED luminosities in this region and eventually they get overcome. 
On the other hand, it should be noticed that in this region the PDF uncertainty for the 
QED-induced channels is almost as large as the luminosity central value, as shown in the 
second inset of the plots in Fig. 9. Thus, by taking into account PDF uncertainties, QED 
luminosities turn out to be compatible at high Mx with QCD luminosities. This feature 
again points to the fact that meaningful comparisons can be pursued only including PDF 
errors, especially when those are large, which is in general true for photon and lepton PDFs. 
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Figure 11: Representative diagrams for partonic processes in e'’'e production. 


Finally, contrary to the luminosities, the 'hjj luminosities maintain the same hierarchy 
all over the range in y. Indeed, the value of y is not directly related to the value of Mx, 
which also in this case is used as factorization scale. Thus, the entire previous argument 
on the suppression of the QCD luminosities with respect to the QED ones does not apply 
for the case of 'hjj. 

In Fig. 9 the value of Mx on the x-axis have been chosen in such a way that similar 
values of r are spanned in the plots for ^/s = 13 TeV and ^/s = 100 TeV. For this reason 
the curves look very similar in the two plots. The only differences are the size of the 68% 
uncertainties, which are larger in the 13 TeV case. This effect is due to the different values 
of the factorisation scale (Mx) for the same values of r, leading to smaller uncertainties 
at 100 TeV, where Mx is larger. As expected, PDF uncertainties are much smaller for 
the QCD initial states (lower inset) than for those involving photons and leptons (central 
inset). 

Similar arguments hold for the comparison of the plot at ^/s = 13 TeV and ^/s = 100 
in Fig. 10. Also in this case the values of Mx^cut, and consequently of Tcut, have been 
chosen in such a way that similar values of r are integrated out in the two cases. With 
higher values of Mx,cut the difference among the luminosities and even their hierarchy can 
change, consistently with what is shown in plots of Fig. 9. However, even for large values 
of Mx,cut, partonic luminosities depend very mildly on y. 

3.2 Production processes 

We now turn to study the effect of lepton PDFs on some key processes at the 13 TeV 
LHC and at the 100 TeV FCC-hh. The analysis of the parton luminosities presented in 
Section 3.1 points to a dependence of the relative size of different parton luminosities on the 
invariant mass of the final-state particles. For this reason, we concentrate on differential 
distributions w.r.t. this kinematic variable. 

All the results have been produced with the help of MadGraph5_aMC@NLO [22] 
using the PDF set C2 in Tab. 1, i.e. apf el_nn23qednlo0118_lept. The relevant SM input 
parameters have been set to the following values: 

as{mz) = 0.118 , Gp = 1.16639 x 10“® , mz = 91.1876 GeV, mw = 80.385 GeV, 
mH = 125GeV, F^ = 2.4952 GeV, F^/= 2.085 GeV. (3.6) 

Since all leptons and quarks, besides the top quark, are considered as partons, they are all 
treated as massless. As a general approach, for all processes we set renormalisation and 


- 15 - 



factorisation scales as fip = ^J-r = Ht/2, where Ht is the scalar sum of the transverse 
masses m'r(*) of the final-state particles, defined as: 

mT{i) = ■ (3.7) 

For the sake of simplicity, all simulations are done at LO and parton level. However, we 
include all the contributions induced by tree-level diagrams, i.e. not only those featuring 
the largest power of ag but also the subleading contributions. Examples will be given for 
the specific processes. 

We remark that here we focus on the impact of the lepton PDFs at the LHC and at 
the FCC-hh. In order to perform a trustworthy similar study for the case of the photon 
PDF, NLO EW corrections cannot be neglected. Notable examples can be found in the 
context of the SM, e.g. in Ref. [23] for the neutral current Drell-Yan or in Ref. [24] for WW 
production, but also in the context of BSM in Ref. [25] for squark-antisquark production. 
Nevertheless, in the following, in order to assess the impact of lepton and photon PDEs, we 
will present results also for photon-induced processes. In general, we will keep separated 
the contributions from initial states with only (anti)quarks and gluons, initial states with 
at least one photon and no leptons, and initial states with at least one lepton. In all the 
plots of this section these three categories are easily identifiable by their distinctive colours, 
which are red, green and blue respectively. 

Neutral Drell-Yan 

We start by considering the case of the production of an electron-positron pair at the LHC 
and at the FCC-hh. Besides the quark-antiquark and photon-photon initial states, this 
process receives, for instance, a new contribution from the electron-positron initial state, 
as can be seen in Fig. 11. Similarly, also and initial states can contribute. 

Fig. 11 suggests that different kinematical distributions can arise from different partonic 
channels such as qq, with only an s-channel diagram, 77 , with t-and n-channel diagrams, 
and e+e“ with s- and t-channel diagrams. It is important to note that all these partonic 
processes yield LO cross sections of thus they all contribute to the same perturbative 

order. 

In Fig. 12 we show the invariant mass distribution of the lepton pair in neutral Drell- 
Yan production at the 13 TeV LHC, with the following cuts: 

PT(e^) > lOGeV, lf?(e^)l < 2.5. (3.8) 

In particular, we consider two invariant mass regions: the low invariant mass region, 
m(e^e“) G [20 GeV, 120 GeV] and the large invariant mass region m(e^e“) > 1 TeV, shown 
in the left and right plot, respectively. In the main panel the qq, 77 and contributions 
are separately displayed with red, green and blue lines, respectively. Similarly, the black 
line is the hadronic cross section, i.e. the sum of the blue, green and red lines. The bands 
displayed around each contribution correspond to the PDF uncertainty at 68 % GL. In both 
the high and low invariant-mass regions we find the contribution of initial-state leptons to 
be fully negligible with respect to the main qq one, with a suppression of more than six and 
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m(e+e-) [GeV] m(e+e-) [GeV] 


Figure 12: e^e production at the LHC. Low invariant mass (left) and high invariant 
mass (right) of the lepton pair. 


more than two orders of magnitude, respectively. As expected, photon-induced processes 
are not negligible for large values of m(e~^e~}. This can be seen also in the first inset of the 
two plots, which display the relative contribution to the hadronic cross section for each of 
the three channels, with £~^£~ too small to be visible in the [1%, 100%] range. In the third 
inset we display the relative size of the 68 % PDF uncertainty band for the three separate 
channels and for the hadronic cross section. The 77 and £~^£~ channels entail much larger 
PDF uncertainties than qq, however the PDF uncertainty of the hadronic cross section 
degrades only at large values of m(£~^£~) due to the 77 contribution. 

In Fig. 13 we show similar plots at the 100 TeV FCC-hh for very large invariant 
masses, m(e~^e~) > 5 TeV. In this case, we also investigate how the results would change 
by loosening the selection cuts down to values that are presumably not achievable at the 
experimental level. In the first row of plots in Fig. 13, we apply only the cuts pT(e^) > 
10 GeV (left plot) and pr(e^) > 100 GeV (right plot). In both plots the contribution 
of the £~^£~ channel is not negligible and in the left plot, for m{e^e~) > 10TeV, it is 
even dominant. This behaviour is due to the e'’'e“ —)■ e'’'e“ partonic cross section with 
massless electrons, which diverges for electrons collinear to the beam pipe. However, once a 
reasonable cut \r]{e^)\ < 4 on the lepton pseudorapidity is set, the contribution of the £'^£~ 
initial state is strongly suppressed. The corresponding distributions with |f?(e^)| < 4 are 
shown in the lower plots of Fig. 13, again with pr(e^) > 10 GeV (left plot) and pr(e^) > 
100 GeV (right plot). In conclusion, although there are kinematic configurations where 
£~^£~ initial states are relevant, they are presumably excluded by standard experimental 
cuts. We observed a similar feature in the pt{^) distributions with m{e^e~) > 2TeV at 
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m(e+e-) [GeV] m(e+e-) [GeV] 




m(e'^e') [GeV] m(e'^e') [GeV] 

Figure 13: e"*“e“ production at the FCC-hh for very high invariant masses with a pT > 
10 GeV cut (left column), a pT > 100 GeV cut (right column), without (top row) and 
with (bottom row) the \r]{e^)\ < 4 cut. 


the LHC; at small pr’s the contribution is dominant, but it is completely suppressed 
if a cut on r]{e^) is imposed. For the 77 initial state, a detailed analysis of the impact of 
the PDF uncertainties at differential level has been performed in Ref. [23]. 
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Dijet 


The next process we consider is dijet production, the dominant process at hadron colliders. 
We look at dijet invariant mass {m{jj)) distributions, for jets with prU) > 10(100) GeV 
at the LHC (FCC-hh) and with or without a cut on the jet pseudo-rapidity {\r]{j)\ < 4) for 
the FCC-hh. At the LHC, instead, we always apply the pseudorapidity cut \r]{j)\ < 2.5. 
We consider two cases: in the first we assume zero probability for leptons and photons 
to fake final state jets, while in second we assume such a probability to be 100%. In 
other words, in the latter case we consider leptons and photons in the final state on the 
same footing of quarks and gluons. The realistic value of the fake rate depends on the 
identity of the particle and on the detector characteristics. At LHC experiments it is 
typically negligible for photons, electrons and muons, while for taus it is about 40%, see 
e.g. Ref. [26]. Moreover, the value of the fake rate is also analysis dependent. Thus, by 
presenting results for the two extreme cases, we provide a conservative estimate of possible 
effects coming from leptons in the initial state. 

Since we consider 2 —)• 2 processes at LO and parton level, the cut on the transverse 
momentum of the final state particles is sufficient to obtain a finite cross section and each 
parton leads to one jet. 

We start by showing plots for the first case (no photons/leptons faking jets), in the 
low and high dijet invariant-mass regions at the LHC and in the very high invariant-mass 
region at the FCC-hh with and without the \r]{j)\ < 4 cut. The LHC and FCC-hh plots are 
shown in Fig. 14 and Fig. 15. In both cases the lepton-initiated contribution is suppressed 
by several orders of magnitude. However, part of the suppression does not originate from 
the PDFs, but from the different perturbative orders of the various partonic channels. The 
gluon-gluon initial state yields 0{ag) terms, which are the dominant ones at LO, while for 
instance the photon-gluon and lepton-lepton initial states yield only 0{asa) and 0 ( 0 ^) 
contributions, respectively. Therefore, the total rates of these subprocesses are suppressed 
also by the subleading order in the perturbative expansion. It is worth mentioning that 
quark-antiquark initial states lead to contributions of C>(a^), O(asa) and which 

we consistently take into account in our analysis. 


The second case, in which we assume that photons and leptons in the final state are 
always reconstructed as jets, is particularly interesting because of new t-channel diagrams 
appearing, e.g., in the initial state. In general, all the £q, I'y or initial states 

give a contribution. The relevant plots for the LHC and the FCC-hh are displayed in 
Fig. 16 and Fig. 17, respectively. Assuming 100% probability for leptons and photons to 
fake final-state jets, the contributions from initial states with at least one lepton (ix) are 
enhanced in such a way that they become as important as those involving photons and no 
leptons. However, both contributions are suppressed by at least four orders of magnitude 
w.r.t. those initiated by QCD partons. 


- 19 




mOj) [GeV] mOj) [GeV] 

Figure 14: Dijet production at the LHC. Low (left) and high (right) invariant mass of 
the dijet pair. Zero probability of photons/leptons faking jets is assumed. 




mOi) [GeV] mOi) [GeV] 

Figure 15: Dijet production at the FCC-hh in the very-high invariant mass region, 
without (left) or with (right) the \ri{j)\ < 4 cut. Zero probability of photons/leptons faking 
jets is assumed. 
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Figure 16: Dijet production at the LHC in the low (left) and high(right) invariant mass 
region. 100% probability of photons/leptons faking jets is assumed. 




mOi) [GeV] mOj) [GeV] 

Figure 17: Dijet production at the FCC-hh in the very-high invariant mass region, 
without (left) or with (right) the \'i]{j)\ < 4 cut. 100% probability of photons/leptons 
faking jets is assumed. 
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W+W and ZZ into 4 leptons 

Next we turn to diboson production at hadron colliders. Firstly, we look at W~^W~ 
production with undecayed W bosons. Then we analyse four-lepton production which is 
the signature emerging from ZZ production with both Z bosons decaying into two leptons. 

In the case of W~^W~ production, we keep the W boson stable so that we can gauge 
the effects due only to the £~^£~ luminosity and not to matrix-element enhancements. The 
latter effects will be instead analysed in the case of ZZ production. Moreover, with stable 
W bosons, it is possible to check if any effect from lepton PDFs could be related to the 
anomaly that has been recently found by the ATLAS Collaboration in the production of 
dibosons with 2 TeV invariant mass in the Run-I of the LHC [27]. 

Differential distributions for the invariant mass m{W~^W~) of the W~^W~ pair are 
shown in Fig. 18, for the LHC in the upper plots and for the FCC-hh in the lower plots. 
For both colliders we consider the case without cuts (plots on the left) and with |r?(VF^)| < 
2.5(4) at the LHC (FCC-hh) (plots on the right). Also for this process the lepton PDF 
contribution can be in general safely neglected. On the contrary, the photon-induced 
processes have to be taken into account for a reliable estimate of the central value of 
the cross section and the associated uncertainty due to PDFs in the large invariant-mass 
region. A similar study at the LHC, including NLO QCD and EW corrections but with no 
PDF uncertainty, have been pursued also in Ref. [24]. The increase of the relative size of 
the 77 -channel for high m{W^W~) is consistent with the behaviour of the and ^qq 
luminosities shown in Fig. 9. In addition, no suppression from s-channel diagrams is present 
in 77 —>■ W'^W~ production, leading to a further relative enhancement with respect to 
the ( 75 -channel at high m{W^W~). On the other hand, in the 77 -channel the W bosons 
are produced more peripherally than in the 55 -channel. Thus, the cut in pseudorapidity 
reduces the relative impact of the 77 channel, but it does not dramatically change the 
qualitative picture. 

In the case of ZZ production, we analyse the four-lepton final state emerging from 
the decay of the Z bosons. Specifically, we consider either four leptons of the same flavour 
(e'''e“e'''e“), or two lepton pairs of different flavours {e~^e~ All diagrams, with and 

without intermediate Z bosons, are taken into account. Thus, all the resonant and non¬ 
resonant contributions have been properly included. It is worth noting that, if we were 
requiring stable Z bosons, no 77 -initiated contribution would appear, since the process 
77 —)• ZZ has no tree-level diagrams. Moreover, as in the case of neutral Drell-Yan, 
final-state leptons collinear to the beam axis lead to divergent cross sections for the £~^£~ 
channel. 

In Figs. 19 and 20 we show differential distributions for the invariant mass of the four- 
lepton system at the LHC and FCC-hh, respectively. In our simulations at the LHC we 
required four leptons with: 


PTi£) > lOGeV, 

\r]{£)\<2.5, 

(3.9) 

while at the FCC-hh they must satisfy: 

Pt{£) > 100 GeV, 

\r]{£)\ <4. 

(3.10) 
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Figure 18: W pair production at the LHC (top row) and at the FCC-hh (bottom row), 
both without cuts (left column) or with a rapidity cut (right column). The rapidity cut is 
\v{W^)\ < 2.5(4) at the LHC (FCC-hh). 


Furthermore, we require that all opposite-sign lepton pairs have an invariant mass > 

20 GeV to avoid collinear singularities due to photon splittings. Again, we clearly observe 
that the .^"’".f^-initiated contribution is completely negligible. The largest relative contribu¬ 
tion to the cross section is at most at the per-mil level at very large values of the four-lepton 
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m(«) [GeV] m(4e) [GeV] 

Figure 19: Four lepton production at the LHC: the case of all four leptons of the same 
flavour (left) and of two different flavour pairs (right) is considered. 




m(«) [GeV] m(4e) [GeV] 

Figure 20: Four lepton production at the FCC-hh: the case of all four leptons of the 
same flavour (left) and of two different flavour pairs (right) is considered. 


invariant mass and tends to become even smaller in the lower invariant-mass region. 
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Same-Sign and/or Different-Flavour lepton pair 

We conclude our review of SM processes discussing two cases which are relevant for both 
SM physics and BSM searches: the production of a same-sign lepton pair and of a different- 
flavour lepton pair. For such searches in the BSM context, the SM background processes 
that are typically considered are those featuring two vector bosons in the final state {ZZ, 
WZ, tt —>• W^W^bb, ...), where the same-sign/different-flavour pair 

of leptons comes together with jets, other (undetected) leptons and/or missing transverse 
energy. Therefore it is worth estimating the lepton-initiated contribution to the SM back¬ 
ground for these processes. 

Before presenting our results, we want to stress one important fact to be kept in mind: 
the introduction of lepton PDFs does not open new production mechanisms. Rather, it 
has to be considered as a different scheme to tackle a given computation, on the same 
footing as the four- and five-flavour schemes in processes involving bottom quarks. For 
example, the process e+e"'' —)■ e+e"'', computed with lepton PDFs, corresponds to the 
process 77 —)• e^e“e'’'e“ where leptons are considered as massive and the e~ degrees of 
freedom are integrated out. This latter scheme has several drawbacks if compared to the 
one with leptons in the initial state: first, the process has a larger final-state multiplicity 
and is therefore more challenging from the computational point of view. On top of this, 
lepton masses are very small as compared to the typical scales probed at the LHC and 
FCC-hh. As a consequence, the numerical evaluation and integration of matrix elements 
can be hampered by instabilities. In particular, large logarithms of the form alog{Q/mi), 
where Q is a typical hard scale of the hadronic process, could spoil the convergence of the 
numerical integration. On the contrary, using lepton PDFs, such logarithms are analytically 
taken into account from m£ up to Qo by the initial condition in Eq. (2.6) and consistently 
resummed from Qo to Q by DGLAP evolution. Therefore, the usage of a scheme with 
lepton PDFs simplifies the computation, and consequently the phenomenological studies, 
for the classes of processes considered here. 

After this clarification, we start by considering the same-sign lepton-pair production. 
We compare the .^"’■^^-initiated case to the process qq —)• W^W^qq, where both the W 
bosons decay leptonically. For this process, which has also been proposed as a way to 
search for and calibrate double-parton scattering at the LHC [28], we do not need to 
consider additional untagged leptons and the cross section is finite even without cuts on 
jets. 

We quote total rates for qq —> W~^W^qq —)■ ViUiqq and for > £~^£^ in Tab. 2. 

For the former, we include all coupling combinations entering at LO, i.e. we evaluate this 
process at the orders \asci^ + ^ As usual, we consider the case of the LHC and of 

the FCC-hh, requiring the two charged leptons to satisfy the cuts of Eqs. (3.9) and (3.10), 
respectively. We can see that in this case the lepton-initiated contribution is not entirely 
negligible as compared to the same-sign W pair production rate, at least at the LHC. In 
this case, the former contribution is 7% of the latter, while this ratio reduces to 1% at the 

^ In Ref. [28] it is shown that the electroweak contribution to the cross section is of the same order of 
magnitude of the QCD one. Thus, it cannot be ignored. 
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£ = e, p,T 

LHG, CT [pb] 

FGG-hh, cT[pb] 

qq —)■ 

3.54 • 10-2 ± 0.4% 
2.64 • 10-3 ± 37% 

6.33 • 10-2 ± 0.3% 
5.86 • 10-^ ± 22% 


Table 2: Rates for same-sign lepton pair production at the LHC and FCC-hh. Leptons 
are required to satisfy the cuts pT > 10(100) GeV, |r/| < 2.5(4) at the LHC (FCC-hh). 
Quoted uncertainties are from PDFs, at 68% CL. 



LHG, cr[pb] 

FGG-hh, cj[pb] 

e'^ IJL~ —>■ 

5.35 • 10-^ ± 38% 

9.93 • 10-3 ± 22% 

e'^T~ —>■ e'^T~ 

5.06 • 10-^ ± 37% 

9.65 • 10-3 ± 22% 

—)• r+/t- 

4.95 • 10-^ ± 38% 

9.46 • 10-3 ± 23% 

pp e^pL~Ue^'^l 

4.90 • 10-3 ± 0.6% 

2.01 • 10-3 ±0.4% 

pp —>■ e~^tr > 20 GeV 

9.96 • 10-2 ± 0.6% 

4.27- 10-2 ±0.5% 


Table 3: Rates for different-flavour lepton-pair production at the LHC and FCC-hh. 
Final-state leptons are required to satisfy the cuts pT > 10(100) GeV, \r]\ < 2.5(4) at the 
LHC (FCC-hh). Quoted uncertainties are from PDFs, at 68% CL. 


FCC-hh. In both cases PDF uncertainties affecting the lepton-initiated contribution are 
quite large. The quoted results are shown for the case of positively charged leptons. In the 
case of negatively charged leptons, the lepton-initiated cross section does not change, while 
the cross section for qq —)• W~W~qq is reduced by approximately a factor 2 (1.4) with 
respect to the positively charged case at the LHC (FCC-hh). The quoted numbers include 
the contributions from all the three leptonic flavours. For qq —)• W^W^qq, different- 
flavour and same-flavour leptons yield two thirds and one third of the total cross section, 
respectively, due to lepton universality. Instead, for —)> the relative contributions 

are 60% and 40% because of non-negligible interference effects that appear only for same- 
flavour and same-sign lepton pair production. 

As far as the case of opposite-sign and different-flavour lepton-pair production is con¬ 
cerned, the total rates for the various flavour combinations are shown in Tab. 3. In the 
table we show lepton-initiated processes and compare them with a non-lepton-initiated 
reaction with the same signature, pp > Leptons are required to satisfy the same 

cuts as in the same-sign case, i.e. Eqs. (3.9) and (3.10). For what concerns pp > 
besides the cuts on leptons we also study the effect of a cut on the missing transverse 
energy $rp < 20 GeV, in order to enhance configuratons with the two leptons back-to-back. 
This cut can be considered realistic for the LHG^, while it may be optimistic at the FGG, 
for which the precise experimental setup is presently unkown. In all cases, lepton-initiated 
contributions are suppressed with respect to pp > and remain below the 1% level 

even with the cut. 


® The typical missing transverse energy resolution in events with a low-pT lepton pair at LHC experiments 
is 10 GeVor less [29-33] 
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Slepton pair it 

The last case that we analyse is relevant for BSM; we consider the hadroproduction of a 
slepton pair (££*) within the Minimal Supersymmetric Standard Model (MSSM). Again, we 
are mainly interested in assessing whether contributions initiated by leptons can affect the 
discovery searches of these particles. Anyway, as done in the previous cases, we consider 
also photon-initiated processes. A similar and more detailed study for 77 initial state 
has been performed in Ref. [25] for the case of squark-antisquark production, where the 77 
contribution has been found to be non-negligible. Here, for slepton-antislepton production, 
we just provide results for total rates. However, at variance with the squark-antisquark 
case, in slepton-antislepton production 77 - as well as -initiated processes are of the 
same perturbative order of qq^ i.e. 0(0^). Therefore, we expect a larger relative impact 
from photon-initiated processes. 

In it production, tree-level Feynman diagrams are very similar to the case of e+e“ 
production. The qq initial state proceeds via (non-resonant) Drell-Yan production, and 
the 77 initial state presents t— and tt—channel diagrams. The i^i~ initial state proceeds 
via both Drell-Yan production and t-channel neutralino exchange. Also, at variance with 
the qq and 77 initial states, at LO the i^i~ initial state can produce iii^ where both the 
chirality and the flavour of ii and £2 different. 

In order to minimise the dependence on the MSSM parameters, we have chosen to 
consider only the cases i = e, n- In this way we avoid the dependence on the left-right 
mixing parameters in the stau sector. Moreover, all the particles besides the selectrons, 
smuons and the lightest neutralino have been decoupled by setting their masses to 5 TeV. 
The neutralino mixing matrix and the physical masses have been calculated at the scale of 
the slepton mass m^, with the help of Suspect [34]. 

In Tab. 4 we list the results for the total cross sections of it production at the LHC 
and FCC-hh, where the soft masses for right-handed and left-handed leptons of the hrst two 
families have been set to m^=200, 500, 1200 GeV and the mass of the lightest neutralino to 
m^o=100 GeV. Physical masses are slightly different from soft masses, however, the small 
differences are completely irrelevant for the study pursued here, so we do not report them 
in the text. The numbers in Tab. 4 refer to the sum of all the possible chirality and flavour 
conhgurations. Also for this process, the contribution due to leptons in the initial state 
is in general negligible. Even setting m^o= 1 GeV, the lepton-lepton contribution results 
slightly enhanced, but the qualitative picture does not change. Gonversely, photon-induced 
processes are relevant and cannot be neglected at the LHG, especially for large masses. 

4 Conclusions and outlook 

In this paper we have presented the first-ever estimate of the lepton content of the proton. 
By implementing the complete LO QED corrections to the DGLAP evolution equations 
in the APFEL program and by means of a model for the lepton PDFs at the initial scale, 
we have produced PDF sets in the LHAPDF6 format containing, besides the usual PDFs 
for quarks and gluons and possibly photons, distributions for all charged leptons. The 
main result of this work are the sets produced according to the ansatz in Eq. (2.6), i.e. 
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LHC, a{pp — 

> £t)[ph] 



Initial state 

mi = 200 GeV 

mi = 500 GeV 

mi = 1200 GeV 

Total 

5.24 • 10-2 ± 4% 

1.31 • 10-3 ± 10% 

1.19-10-3 ±38% 

m 

4.98 • 10-2 ± 1% 

1.18- 10-3 ±2% 

7.52 • 10-® ± 6% 

77 

2.62 • 10-3 ± 72% 

1.32-10-^ ±84% 

4.38 • 10-® ± 95% 


6.70 • lO-’^ ± 92% 

4.92 • 10-3 ± 97% 

1.68-10-9 ±99% 


FCC-hh, a{pp rt)[pb] 


Initial state 

mi = 200 GeV 

mi = 500 GeV 

mi = 1200 GeV 

Total 

8.42 • 10-3 ± 2% 

4.88 • 10-2 ± 3% 

2.35 • 10-3 ± 5% 

qq 

8.13 • 10-3 ± 1% 

4.64 • 10-2 ± 1% 

2.19 • 10-3 ± 1% 

77 

2.94 • 10-2 ± 46% 

2.32 • 10-3 ± 52% 

1.58-10-^ ±62% 

e+i- 

6.22 • 10-® ± 74% 

8.09 • 10-^ ± 80% 

7.86 • 10-3 ± 87% 


Table 4: Total rates for slepton-pair production at the LHC and at the FCC-hh. Quoted 
uncertainties refer to PDF uncertainties at the 68 % CL. 


by assuming that leptons are produced in pairs by photon splitting at the respective mass 
scales. 

By studying the feature of the PDFs, we found that our results are consistent with the 
expectations. Lepton PDFs are strongly correlated to the photon PDF and suppressed by 
a factor a w.r.t. it. As a consequence of their small size, their phenomenological impact 
at hadron colliders is likely to be limited. Also, the uncertainties of lepton PDFs can be 
very large, as in the case of the photon PDF. We explicitly checked the phenomenological 
impact of lepton-induced partonic channels for several relevant SM processes at the LHC 
at 13 TeV and at the 100 TeV FCC-hh. With reasonable experimental cuts, the contribu¬ 
tions of lepton PDFs is usually extremely suppressed. On the contrary, the photon PDF 
contribution is not negligible and has to be carefully studied for precise predictions at the 
LHC and FCC-hh. We also presented representative results for BSM processes (slepton- 
antislepton production in the MSSM) and rare SM processes such as same-sign and/or 
different-flavour dilepton production. More detailed studies, including higher-order effects 
as well as realistic cuts and background simulations, are envisaged, since the role of the 
lepton PDFs could be relevant in this context. 

The combined QCD-I-QED evolution with lepton PDFs has been implemented in APFEL 
version 2.4.0 and later, publicly available at: 

https://github.com/scarrazza/apfel 

The sets of PDFs used in this work, produced with APFEL, are also publicly available 
in the LHAPDF6 library format from the APFEL webpage^: 

apf el Jiin23nlo0118_lep0, apf el Jnn23niilo0118_lep0, 

^http://apfel.hepforge.org/ 
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apfel_nn23qedlo0118_lept0, apfel_nn23qednlo0118_lept0, apfel_nn23qednnlo0118_lept0, 
apfel_nn23qedlo0118_lept, apfel_nn23qednlo0118_lept, apfeljnn23qednnlo0118_lept, 
apfel_mrst04qed_lept0, apfel_mrst04qed_lept 
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A The combined QCD+QED evolution in the presence of lepton PDFs 

The DGLAP evolution equations employed in this work for the unified QCDj-QED evo¬ 
lution in the presence of lepton PDFs consist of a system of coupled differential equations 
which satisfy both QCD and QED evolution equations. In order to simultaneously diago- 
nalise as much as possible the respective evolution matrices, avoiding unnecessary couplings 
between parton distributions, we have adopted the following evolution basis: 

1 : 9 

2 : 7 

3 : S = Yiu + Tifi 

4 : As = - Srf 

5 : 

6 : Tf = ii+ - c+ 

7 : r 2 “ = ii+ c+ - 2t+ 

8 : Tf = d+- s+ 

9 : = d-^ + s-^ - 26+ 

10 : Ti = If - 1+ 

11 : Ti = lt + l+-2lt 

where we have defined: 

g+ = g ± g with q = u,d, s, c, b, t, 

(A.2) 

lf=i~± £+ with i = e, p,, r. 


12 

13 

14 

15 

16 

17 

18 

19 

20 


V = Vu + Vd 
Ay = Vu-Vd 
Ve 

= u~ — c~ 

Vf- = u~ + c~ — 2t~ 
Vf = d- -s- 

= d- +s- - 2b- 

vi = if-i-, 

= 17 +17- 2K 


(AT) 
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and: 


5^11 — + c~*~ + i~*~, Vu — u ~\~ c + t , 

Yjd = d~*~ + s~*~ + 6~*~, Vd = d + s -\- b , (-^-3) 

= /+ + /+ + Vi = . 

When considering LO QED corrections to the DGLAP equations, the evolution equa¬ 
tions in the basis given in Eq. (A.l) decouple into different sectors that evolve independently 
as shown below. It should also be noticed that considering only LO QED corrections al¬ 
lows one to completely separate QCD and QED corrections. Indeed, defining P the QCD 
splitting functions (up to any of the known orders in q;<j [35-39]) and P = {a/47r)P^^'^ the 
LO QED splitting functions, with P^^'^ given in Eq. (2.3), we have: 


• the singlet sector: 


/ g\ 



/ 

p 

^99 

0 

p 

^99 

0 o\ 

7 




0 

0 

0 

0 0 

S 

= 



‘^'’^fPq9 

0 

p 

^qq 

0 0 

As 



n 


0 

riu—nd i D 
n/ Wot 

- P+) P+ 0 

vsj 



V 

0 

0 

0 

0 oy 


/O 0 0 0 0 \ 




0 7]^ P-yq Tj P'^q P'jq 



7 

0 2 e|.P 5 .^ V'^Pqq 'n~Pqq 0 


(8) 

s 

0 2dlPqy g-pqq ri+Pqq 0 



As 

\0 2nePqy 0 0 PqgJ 





(A.4) 


• the quark valence sector: 




2 


d 

dy? 



riu-nd 


P^ 0 \ 

(pV _ p-) p- I 


/ ri+Pgq IJ Pqg\ 

\V~Pqq V^Pqq) 



(A.5) 


• the lepton valence sector: 


ndVi 


(A.6) 
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• the quark non-singlet sector: 


= iP^ + elP,,)0Tr^2, 


dp? 


„ dT?r, _ , 

= (J’+ + e3J’„)®r^2. 


?' Q .,2 




dp? 

Mtp 

dp^ 


= {P-+elP,,)^Vr,2, 


= {P-+ejP,g)®Vf^,, 


(A.7) 


and the lepton non-singlet sector; 


dT^ 

,,2 3,8 _ p p. rpl 

^ dp^ ~-^ 3,8 5 


2 ^^3^8 - / 

p3 = Pm®yifi^ 


dp^ 


(A.8) 


where and are the electric charges of the up- and down-type quarks, nf is the number 
of of active quark flavours, and are the number of up- and down-type active quark 
flavours (such that n/ = Uu + Ud) and ni is the number of active leptons. Moreover, we 
have defined: 

e| = Nciriuel + ndej ), 

= Nciuuel - Udej) , (A.9) 

^ {el ± ej) , 

being Nc the number of colours. 

The basis presented here is implemented in APFEL and accessible through the QUniD 
solution. 
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